A late-accelerating universe with no dark energy 
— and a finite-temperature big bang 
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Brane-world models offer the possibility of explaining the late acceleration of the universe via 
infra-red modifications to General Relativity, rather than a dark energy field. However, one also 
expects ultra-violet modifications to General Relativity, when high-energy stringy effects in the early 
universe begin to grow. We generalize the DGP brane-world model via an ultra-violet modification, 
in the form of a Gauss-Bonnet term in the bulk action. The combination of infra-red and ultra-violet 
modifications produces an intriguing cosmology. The DGP feature of late-time acceleration without 
dark energy is preserved, but there is an entirely new feature - there is no infinite-temperature big 
bang in the early universe. The universe starts with finite density and pressure, from a "sudden" 
curvature singularity. 



INTRODUCTION 



The standard cosmology based on General Relativ- 
ity and inflation has been remarkably successful. But 
there remain deep puzzles left for theorists to resolve 
- what is the cause of the late-time acceleration of the 
universe (the "dark energy" problem)? how is the clas- 
sical big bang singularity removed by quantum gravity 
effects? One approach to start tackling these problems 
is via the brane-world scenario, which is motivated by 
string theory. Most brane-world models, including those 
of Randall-Sundrum type Q , produce ultra-violet mod- 
ifications to General Relativity, with extra-dimensional 
gravity dominating at high energies. However it is also 
possible for extra-dimensional gravity to dominate at low 
energies, leading to infra-red modifications of General 
RelativiW 0,13. The Dvali-Gabadadze-Porrati (DGP) 
models Q (see also 0) achieve this via a brane induced 
gravity effect. 

The generalization of the DGP models to cosmology 
lead to late-accelerating cosmologies, even in the absence 
of a dark energy field This exciting feature of "self- 
acceleration" may help towards a new resolution to the 
dark energy problem. But the models suffer from the 
short-coming that they do not modify 4D gravity at high 
energies, where we expect stringy corrections to start 
having an effect. How can we generalize the DGP models 
so that they also show ultra-violet modifications to Gen- 
eral Relativity? One possibility is to introduce into the 
gravitational action a term that is associated with higher- 
energy stringy corrections - the Gauss-Bonnet term 
Indeed, in certain realizations of string theory, the ghost- 
free GB term in the bulk action may naturally lead to a 
DGP induced gravity term on the brane boundary 0. 

We investigate what happens when a Gauss-Bonnet 
term is introduced in the 5D Minkowski bulk containing 
a Friedmann brane with DGP induced gravity. As we 
will show, this combination of infra-red and ultra-violet 
modifications leads to intriguing cosmological models. 



2. FIELD EQUATIONS 

The gravitational action contains the Gauss-Bonnet 
(GB) term in the bulk, as a correction to the Einstcin- 
Hilbert term, and the Induced Gravity (IG) term on the 
brane, 

+ / d^xy^i?^, (1) 

^'^S J brane 

where a (> 0) is the GB coupling constant [8| and r (> 0) 
is the IG "cross-over" scale, which marks the transition 
from 4D to 5D gravity. The DGP models are the special 
case a = 0, and in this case the cross-over scale defines 
an effective 4D gravitational constant via k1 = K^/r . 

We assume mirror {Z2) symmetry about the brane. 
The standard energy conservation equation holds on the 
brane. 



p + 3H{l + w)p = 0, 



P/P- 



(2) 



The modified Friedmann equation was found in the most 
general case (where the bulk contains a black hole and 
a cosmological constant, and the brane has tension) in 
Rcf. 0. For a spatially flat brane without tension, in a 
Minkowski bulk, 



l-f -a ( i/^ + - 
3 V 2 



(i?^ - $) = 



where $ is determined by 

$ -I- 2a$2 ^ . 



,(3) 



(4) 



(In the most general case, the right-hand side of this 
condition is nonzero 0.) Equation Q has solutions $ = 
0, — l/2a, but here we only consider $ = 0, since the 
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FIG. 1: DGP and GB solutions of the Friedmann equation 
{H vs p) for a Minkowski bulk. (Arrows indicate the fiow of 
brane proper time t, with t = oo at p = 0.) 



second solution has no IG limit and thus docs not include 
the DGP model [l^. 

The DGP models have a = <f> = and the Friedmann 
equation (j^J reduces to a quadratic in with solutions 



3. DGP BRANE WITH GB BULK GRAVITY: 
COMBINING UV AND IR MODIFICATIONS 

The DGP(+) models are attractive for cosmology since 
they accelerate at late times, without the need for dark 
energy, when gravity begins to leak off the brane, i.e., 
when the 5D Ricci term in Eq. begins to dominate 
over the 4D Ricci term. At early times, the 4D term 
dominates and General Relativity is recovered (in the 
background - note that the perturbations are not Gen- 
eral Relativistic ^3). The DGP models are in some sense 
"unbalanced" , since they do not include ultra-violet mod- 
ifications to cosmological dynamics. In order to modify 
4D gravity at high energies as well as low energies, we 
can include a GB term in the action. 

The combined Gauss-Bonnet Induced Gravity (GBIG) 
model has a DGP brane in a Minkowski bulk with 
Einstein-Gauss-Bonnet gravity. The GBIG Friedmann 
equation follows from putting $ = in Eq. ©. Defining 
dimcnsionlcss variables, 



8a , rui t 

7 = TT' h = Hr, ^ = —p , T = - 
6r^ 6 r 

the GBIG Friedmann equation becomes 



(12) 



(13) 



2 ni 
±-H+^p. 
r 6r 



(5) 



There are two branches DGP(±) for the two signs on the 
right (corresponding to different embeddings of the brane 
in the Minkowski bulk). Both branches have a 4D limit 
at high energies, 

DGP(±): H:$>r-^ ^ (x p , (6) 

while at low energies. 



DGP(+) 
DGP(-) 
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P 



r 
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(7) 
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DGP(-) has a non-standard (and non-accelerating) late 
universe. The self-accelerating DGP(4-) branch is of most 
interest for cosmology, and we focus here on this model 
and its generalization via a GB term. 

The pure GB model with a Minkowski bulk has r = 
$ = and the Friedmann equation Q reduces to 



1 



5 „2 



(9) 



which is a cubic in H^. This GB Friedmann equation has 
no 4D limit: 



GB high energy: 
GB low energy: 



H » a- 
H « a- 



cx p2/3 , (10) 
(x p^ . (11) 



The Friedmann equations for pure DGP and pure GB 
models with a Minkowski bulk are compared in Fig. ^ 



while the conservation equation becomes 
p.' + 3h(l + w)p = 0, 



(14) 



where a dash denotes d/dr, and h ~ a' /a. 

Combining Eqs. H13() and ()14|l . we find the modified 
Raychaudhuri equation. 



h' 



3/i(l + w){h?' — fi) 



A{-fh^ + l)(37/i2 + 1) - 2(/72 - p) ■ 
The acceleration a" /a = h' + h'^ is then given by 



(15) 



_ Ah\jh'^ + l)(37/i2 + 1) - (fe2 _ ^)[2fe2 - 3(1 + w)p] 
a ~ 4{-fh^ + l){3-fh^ + 1) - 2{h^ - p) 

. 

The GB correction, via a non-zero value of 7, intro- 
duces significant complexity to the Friedmann equation, 
which becomes cubic in h'^, as opposed to the quadratic 
DGP(±) case, 7 = 0, for which 



P 



2±2^p+l 



(17) 



This additional complexity has a dramatic effect on the 
dynamics of the DGP(+) model, as shown in Fig.|21 The 
contribution of GB gravity at early times removes the 
infinite density big bang, and the universe starts at finite 
maximum density and finite pressure (but, as we show 
below, with infinite curvature). Furthermore, there are 
two such solutions, each with late-time self-acceleration, 
marked GBIGl and 2 on the plots. Since GBIG2 is accel- 
erating throughout its evolution (actually supcr-infiating. 
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DGP(t) 



FIG. 2: Solutions of the Friedmann equation {h vs /i) for the 
DGP(+) model and its Gauss-Bonnet corrections, GBIGl and 
GBIG2. The curves are independent of the equation of state 
w. Arrows indicate the flow of normalized brane proper time 
r, with r = 00 at = 0. (Here 7 = 0.05.) 




big bang, p = 00, in common with the DGP(±) and 
GB models in Fig. ^ By contrast, GBIGl-2 have a 
finite-temperature big bang, since the density is bounded 
above, 



(18) 



where fii (which is positive only for 7 < 1/16), is found 
below, in Eq. (|^ . 

The finite-density beginning was pointed out in 
Ref. 0, where the cubic for the general case (i.e., with 
tension, bulk cosmological constant and bulk black hole) 
was qualitatively analyzed. Here, we focus on the sim- 
plest generalization of the DGP models, and give a de- 
tailed quantitative analysis of the cosmological dynam- 
ics. In particular, our analysis shows that one solution 
is not bounded, which was not noticed in Ref. The 
numerical plots of the Friedmann equation in Figs.|21and 
13 are crucial to a proper understanding of the algebraic 
analysis of the cubic roots. 

A detailed analysis of the cubic equation (|13|) con- 
firms the numerical results, and shows that (for /i > 0): 



< 7 < — : 3 real roots, GBIGl-3, 
16 

7 > — : 1 real root, GBIG3. 
' - 16 

The real roots arc given as follows: 

• For < 7 < 1/16; the roots GBIGl-2 are 



(19) 
(20) 



1 - 87 



4-^2^2 ^ + 2^/^003(6* + — ) for A* < Mi , (21) 
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where n = 4 for GBIGl, n = 2 for GBIG2, and the root 
GBIG3 is 



2u2 1-87 



2\/^Qcos6' for /i < /ii. 



(22) 



FIG. 3: The DGP(-) model and its GB correction, GBIG3. 



S+ + S- for /-t > Mi . 
• For 7 > 1/16: the only real root GBIG3 is 

472/12 = i—^ + S+ + S-. (23) 



h' > 0), the physically relevant self-accelerating solution 
is GBIGl. 

The cubic in h^, Eq. (|13|l . has three real roots when < 
7 < 1/16 (sec below). Two of these roots correspond to 
GBIGl-2, which are modifications of the DGP(-I-) model. 
The third root GBIG3 is a modification of the DGP(-) 
model, as illustrated in Fig. 13 Note that the curves in 
these figures are independent of the equation of state w 
of the matter content of the universe - w will determine 
the time evolution of the universe along the curves, via 
the conservation equation (|14|l . 

The plots show that GBIG3 starts with a standard 



In the above, S±, Q, R, 6 arc defined by 

1/3 



S- 



± 



R ± v/i?2 + Q3 
1 



Q=^(M + 2)--(l-87)^ 



(24) 
(25) 



R = 87V^ - - 
cos 36* = i?// 



47^ (1 - 87)^ 

^(l-87)(M + 2) + ^^y^,(26) 



(27) 



The explicit form of the solutions can be used to con- 
firm the features in Figs. [3 and |21 Equations if^ - lf^ 
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show that GBIG3 starts with a big bang, h, fi —f oo, with 
~ /i^/'^ near the big bang. This is the same as the high- 
energy behaviour of the pure GB model, Eq. ((Tn|l - the 
GB effect dominates at high energies in GBIG3. This is 
not the case for GBIGl-2. 

The maximum density feature of GBIGl-2 is more eas- 
ily confirmed by analysing the turning points of /i as a 
function of h^. The Fricdmann equation H13|) gives 



/i-2(7/i2-f 1) (37/12 -^1) 



Substituting dfj./d{h'^) = into Eq. (O, we find that 



(28) 



30- h. 



hi 
Mi 



1±V1-127 
67 ' 

±hf (2VT^T2^T1 



(29) 
(30) 



The second equation shows that positive maximum den- 
sity only arises for the upper sign and with 7 < 1/16, 
in agreement with the cubic analysis. Thus the initial 
Hubble rate and density for GBIGl-2 are 



0.01 0.02 0.03 0.04 0.05 



FIG. 4; The dependence in GBIGl-2 of the initial expansion 
rate and density on 7. 



hi = 



Mi 



1 + yi - 127 
67 ' 

i/^f (2^1^^- 1 



(31) 
(32) 



If 7 = 0, then GBIGl-2 reduce to DGP(+), and 
fii — 00. Note that 



hi>A. 



(33) 



The case 7 = 1/16, /i^ = corresponds to a vacuum brane 
with de Sitter expansion, and h = hi = A, generalizing 
the DGP(-I-) vacuum dc Sitter solution |^. 

The late-time asymptotic value of the expansion rate, 
as jU 0, is 



hr. 



2V2 



7 



1 - 87 T VI-I67 



1/2 



(34) 



where the minus sign corresponds to GBIGl and the plus 
sign to GBIG2. In the hmit 7 -> 0, GBIGl recovers the 
DGP(+) case, h^ = 2, while for GBIG2, 00; 
the parabolic GBIGl-2 curve in Fig. [21 "unwraps" and 
transforms into the DGP(-I-) curve. Equations 119() and 
show that 



2 < /loo < 4 for GBIGl, 



(35) 



while 4 < /loo < 00 for GBIG2. 

The behaviour of the key GBIGl-2 parameters is il- 
lustrated in Figs. 21 and El 
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FIG. 5: The GBIGl-2 late-time asymptotic expansion rate 
as a function of 7. 



holds. By Eq. (|12|l . this means that the GB length scale 
Lgb = \foL must be below a maximum threshold deter- 
mined by the IG cross-over scale: 



^<16 ^^^^ 



r- 1 /3 
I OL < —\ — r . 

8 V 2 



(36) 



4. COSMOLOGICAL DYNAMICS 

The GBIGl model, which is the physically relevant 
generalization of the DGP(+) model, exists if Eq. H19|) 



If the GB term is taken as the correction term in certain 



SO 



string theories, then Lgb ~ -^^string, while r ^.iq 
that this bound is easily satisfied. 

When Eq. H36|) holds, the universe starts with a max- 
imum density pi and maximum Hubble rate Hi, and 
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evolves to an asymptotic vacuum de Sitter state: 




(37) 
(38) 



At any epoch tq, the proper time back to the beginning 
is 



To - n 



da 
ah 



(39) 



Since a and h are nonzero on the interval of integration, 
the time back to the beginning is finite. 

The current Hubble rate can be approximated by the fi- 
nal de Sitter Hubble rate, Hq ~ i?oo, so that by Eq. (|35|l . 
the cross-over scale obeys 



(40) 



In the DGP(-H) hmit, r ~ 2H^^. The effect of GB gravity 
is to increase r but not beyond r ~ 47?^^. 

However, there is a UV-IR "bootstrap" operating to 
severely limit the GB effect at late times. The key point 
is that appreciable late-time GB effects require an in- 
crease in 7, whereas the primordial Hubble rate Hi is 
suppressed by an increase in 7 - as shown in Figs. ^ and 
m Equations ||2U and imply that 



Hi » Ho 



7 < 



f 

16 



(41) 



Thus the GBIGl model does not alleviate the DGP(^ 
fine-tuning problem of a very large cross-over scale, r 
H^^ - (10-33 eV)-i. 

The GBIGl Friedman equation H21|l gives 

1 - 87 



H^ = 



I27V2 

V87^(rK^p + 6)-(l~87)^ 



, , 47r 



(42) 



where 



cos 30 (p) = 

2167-^At^ - 367^(1 - 87)(/^ -f 2) + (1 - 87)3 
[(l-87)2-2472(/x-h2)]3/2 



(43) 



A more convenient form of the Fricdmann equation fol- 
lows from solving Eq. H13|l for /i, 



fi = - 2h{"fh^ + 1) , hoo <h <hi. 



(44) 



By expanding to first order in — h?^, we find that at 
late times, 



(45) 



Taking the DGP(4-) limit hoc ^ 2, and comparing with 
Eqs. lO and 117|) . we find that the effective Newton con- 
stant in GBIGl is 



(46) 



where G5 = KrjSir is the fundamental, 5D gravitational 
constant. In the DGP(-I-) case G = G^/r. 

Equation H46|l gives a relation for the fundamental 
Planck scale A/5 



Ml 



rHn 



rHn 



p 



(47) 



where Mp is the effective 4D Planck scale, and we used 
Hoo ^ Hq. As r 4i7g"^ (its upper limit), so M5 in- 
creases. This is very different from the DGP(-I-) case, 
where M| = Afp/r, so that is constrained to be very 
low, < 100 MeV. In principle, GB gravity allows us 
to solve the problem of a very low fundamental Planck 
scale - but in practice the UV-IR bootstrap, Eq. H41|l . 
means that 7 ~ so that is effectively the same as 
in the DGP(-I-) case. 

What is the nature of the beginning of the universe in 
GBIGl? We can use Eq. ^ in Eq. for matter with 
w > —1, to analyze the initial state, d^/d{h'^) 0+. We 
find that h[ = —00, i.e., infinite deceleration, 

< = -00 . (48) 

(For GBIG2, with dfi/d{h'^) 0-, we have h[ = +00.) 
The initial state has no infinite-temperature big bang, 
but it has infinite deceleration, and thus infinite Ricci 
curvature. The brane universe is born in a "quiescent" 
singularity. (Although similar singularities may be found 
in Induced Gravity models , they arise from the spe- 
cial extra effect of a bulk black hole or a negative brane 
tension.) The key point is that neither the DGP(-f) 
model nor the GB model avoid the standard big bang, as 
shown in Eqs. © and ((Tn|l . But together, the IG and GB 
effects combine in a "nonlinear" way to produce entirely 
new behaviour. If we switch off either of these effects, 
the infinite-temperature big bang reappears. 

This singularity is reminiscent of the "sudden" (future) 
singularities in General Relativity - but unlike those 
singularities, the GBIGl-2 singularity has finite pressure. 
The initial curvature singularity signals a breakdown of 
the brane spacetime. The (Minkowski) bulk remains reg- 
ular, but the imbedding of the brane becomes singular. 
Higher-order quantum-gravity effects will be needed to 
cure this singularity. The fact that the matter is regular 
at the singularity indicates that the singularity is weaker 
than a standard big bang singularity, and may be easier 
to "cure" with quantum corrections. 

By performing an expansion near the initial state, us- 
ing Eq. H44|) . we find that the primordial Hubble rate in 
GBIGl, after the infinite deceleration at the birth of the 
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have used the effective equation of state 




FIG. 6: The acceleration / = a"/a vs x = , for a GBIGl 
cosmology with inflation, followed by radiation domination, 
followed by late-time self-acceleration. Arrows indicate the 
flow of brane proper time. Here 7 = 0.05, and n — 0.8 in 
Eq. (inj- 



universe, is given by 



3Vr2 - 32a 



1/2 



(Pi - pY^'- (49) 



This is independent of the equation of state w. If there 
is primordial inflation in the GBIGl universe, then the 
acceleration a" will become positive. For a realistic 
model (satisfying nucleosynthesis and other constraints), 
a" must subsequently become negative again, so that the 
universe decelerates during radiation- and early matter- 
domination. Finally, a" will become positive again as the 
late universe self-accelerates. 

We can simplify the expression H16|) for the acceleration 
in GBIGl via Eq. (gH), 



/ = 



2x{3'yx + 1- y/x) + 3(1 + w)[xy/x - 2x{-fx + 1)] 



2(373; + 1 - v^) 

(50) 

where / = a"/a,x = h^. For a given w(x), we can 
plot f{x). We show an example in Fig. |S1 of a simple 
model, with primordial inflation followed by radiation 
domination, followed by late-time self-acceleration. We 



-0.9 n{h? ~ hl^) + hl^ < X < hf , , . 
1/3 hl,<x< 7i{hf - hi,) + hi, . ^'^^ 



Here < n < 1 is a parameter determining the time of 
reheating (with n = corresponding to no inflation and 
n = 1 to no reheating/ radiation). 

5. CONCLUSIONS 

In summary, the GBIGl model provides an intriguing 
generalization of the DGP(+) model - the Gauss-Bonnet 
(ultra-violet) correction to the (infra-red) Induced Grav- 
ity preserves the late-time self-acceleration of the uni- 
verse, but leads to striking new behaviour in the early 
universe. Although there is still a curvature singularity 
at the beginning, the density, pressure and temperature 
arc finite. This model deserves further investigation as a 
viable cosmological model. Future work will impose con- 
straints on the model parameters from nucleosynthesis 
and Supernova redshifts (compare Ref. 0|). We expect 
that these constraints will not differ appreciably from 
the DGP(+) case, given the very small value of the GB 
parameter 7 that is imposed by the UV-IR bootstrap. 
However, a nonzero 7, no matter how small, leads to 
dramatic and nonperturbative changes at high energies 
in the primordial universe. 

From a theoretical viewpoint, it will also be impor- 
tant to investigate how the GB term affects the issues 
of strong coupling and ghosts in the DGP(-I-) cosmologi- 
cal model [13 . Can the GB term provide a lowest-order 
ultra-violet completion of the DGP(-f ) theory? The anal- 
ysis of perturbations about a Minkowski brane with GB 
and IG terms ^3 gives a starting point for tackling the 
Friedmann brane case. 
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